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O . Abstract 

We investigate the interaction between H-dibaryons employing a quark clus- 
ter model with a one-gluon-exchange potential and an effective meson exchange 
potential (EMEP). A deeply-bound state of two H-dibaryons due to the medium 
range attraction of the EMEP is obtained. The bound H-H system has a size 



of about 0.8 ~ 0.9 fm because of the short-range repulsion generated by the 
color-magnetic interaction and the Pauli principle. 

1 Introduction 

The H-dibaryon is a doubly strange six-quark state with spin and isospin 0. Since its 
first prediction in 1977 flT], it has been the subject of many theoretical and experimental 
studies as a possible candidate of a strongly bound exotic state. The main question 
is whether its mass is lighter than the AA threshold, because such a light H-dibaryon 
is stable against strong decay. However, theoretical estimates of its mass have varied 
widely ranging from deeply bound to unbound, depending on the model or theory used 
for the calculations 0. Experimentally, there is no conclusive result yet, although 
several candidates for an H-dibaryon decay have been reported ||. Many efforts to 
search for the H-dibaryon are under way H and planned. A double-hypernucleus 
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event with a AA binding energy less than about 30 MeV, has recently been found in 
an emulsion experiment at KEK J7|. This excludes the possibility of a deeply-bound 
H-dibaryon. 

The main dynamical factor responsible for a deeply-bound H-dibaryon is the color- 
magnetic interaction (CMI) in the one-gluon-exchange potential (OGEP). The CMI 
plays an important role in explaining the low-lying hadron spectra and the short-range 
part of the baryon-baryon interaction. In the H-dibaryon, the CMI is attractive. On 
the other hand, it acts as a repulsive force for the nucleon-nucleon interaction and also 
for almost all nucleon-hyperon and hyperon-hyperon interactions. As is shown later, 
the CMI also generates a repulsive force between two H-dibaryons. Another important 
factor determining the quark dynamics is the Pauli exclusion principle. It is known 
that the Pauli principle gives repulsion for some channels of the nucleon-hyperon and 
hyperon-hyperon systems. Furthermore, because of the color exchange character of 
the OGEP, the CMI is only effective when the quark exchange interaction between 
different baryons is taken into account. 

Quark exchanges between clusters can be correctly treated in the quark cluster 
model through the antisymmetrization of the wave function. Therefore, we employ the 
quark cluster model to study the interaction between H-dibaryons in this paper. 

Assuming that the CMI plays a key role in determining the properties of the H-H 
interaction, Tamagaki suggested the possibility that a new form of hadronic matter, 
called H-matter, appears at densities several times higher than normal nuclear density 
II. In this pioneering study, a simple H-H interaction model consisting of a hard-core 
potential plus a square well attractive potential outside the core was used, because no 
microscopic calculation of the H-H interaction was available at that time. 

In the present study, we extract more quantitative information on the H-H inter- 
action by employing the quark cluster model used in Ref. ||. In this model, quarks are 
confined by an inter-quark potential and subject to the OGEP |J. The medium range 
attraction is taken into account by an effective meson exchange potential (EMEP) be- 
tween H-dibaryons. Here we take only flavour singlet scalar meson exchange between 
the H-dibaryons into account. In the next section, we briefly explain the model used 
in the calculation and give the necessary matrix elements of the spin-flavor-color op- 
erators. The results of the calculation and the discussion of the properties of the H-H 
interaction are given in section |[ Section f| is devoted to a summary. 
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2 Model 



In the present quark cluster model, the H-dibaryon wave function (/)($) is assumed to 
be a (Os) 6 configuration divided by the Os wave function of the center-of-mass of the 



H-dibaryon with the size parameter y 1/66 for the orbital part, where £ is the 

internal coordinate. 

<K£) = ^)«S([222] SF )C([222] C ). (1) 

Here, <S([222]sf) and C([222]c) are the spin-flavor and color parts, respectively, with 
their permutational symmetries in parentheses. The total wave function of the H H 
system is written as 

^(^^2,r)=Am i m 2 )x(r)], (2) 
where the wave function of the relative motion between the H-dibaryons x{ r ) is intro- 
duced and r is the relative coordinate between the H-dibaryons. The wave function of 
the relative motion is further decomposed into components of the angular momentum 
of the relative motion, L, as 



x(r) = £*V^W*)- (3) 

LM r 

In this paper, only relative S-waves are considered. 

The antisymmetrization operator A can be expressed as 

A=l- 36P 67 + 225P 57 P 68 - 400P 47 P 58 P 69 , (4) 

where is the permutation operator. Because each H-dibaryon consists of 6 quarks, 
there are 0,1,2 and 3 quark exchanges between H-dibaryons which are shown in Fig. 
[I]. Because we symmetrize the two H-dibaryons system under a permutation of the 
two H-dibaryons, 4,5 and 6 quark exchanges are taken into account by the 0,1,2 and 3 
quark exchanges and the exchange of the clusters as a whole. 

The Hamiltonian H consists of the kinetic energy term, K, and the interaction 
term, V. 

H = K + V, (5) 

i=i 

where Kq is the center-of-mass energy of the twelve-quark system, and 

^ = E(^ OGEP + ^ onf )- (7) 

i<j 



Here, V® GEP and V™ are the one-gluon-exchange potential (OGEP) and the confine- 
ment potential, respectively: 



yOGEP _ ^£ 

ij A 



1 7T N A 2 

Ofaj) 1 + -<T; • cr- 

\_Tij rriirrij \ 6 



\ ■ Aj, (8) 



V^ onf = -acrJ-Ai • A r (9) 

The term containing er^ • cr^Aj • Xj in eq.(^j) is known as the color-magnetic interac- 
tion (CMI) term. It is important for the short-range behavior of the baryon-baryon 
interaction and also plays a crucial role in the short-range H H interaction. 

Flavor SU(3) symmetry breaking in the CMI is taken into account by using the 
parametrization: 

7T /rriirrij -> fan /ml, 

where fa = 1 when both i and j are not an s-quark, fa = £i when either i or j is 
an s-quark, and fa = £2 when both % and j are s-quarks. Note that we take £2 as an 
independent parameter, instead of being £f, which is straightforwardly expected from 
eq.(||). As in ref.||, we use here the parameter set which reproduces the experimental 
values of the mass splittings among AA, NS and EE correctly. 

The other parameters entering OGEP are determined from the experimental masses 
of octet and decuplet baryons || . Their values are reproduced here in Table [I]. 

The equation that determines the relative wave function x(r) is called the Resonat- 
ing Group Method (RGM) equation. 

/ {Hrgm(t, r') - EN RGM (r, r')} X (r')dr' = 0, (10) 

where A^rgm (-Hrgm) is the normalization and Hamiltonian kernel. These kernels are 
calculated as matrix elements of 1 and H using the wave function of eq.(^). 

Details on how to evaluate the normalization and Hamiltonian kernels can be found, 
for example, in refs. [pi], |TT |. Here, we show the spin-flavor-color (SFC) parts of the 



matrix elements. Because the SFC parts of the H-dibaryon wave functions are totally 
antisymmetrized, it is sufficient to calculate the following matrix elements: 

(IfM/fM j e .. A . . Ai(or £.. A . . Xj(7 . . a^pjr™) | #(™)#(™)) 



{HWHW I iijXi ■ Xj(ov ■ ■ <T 3 )Pi° Tc) Pil Tc) Pt c) I H^H^). (11) 

Here, | H (aTC ^H (aTC ^) is the SFC part of the shell model wave function in eq. (§) . There 
are 66 pairs interacting via the two-body interaction, V^-, and they are classified into 
ten topologically different types of graphs as illustrated in Fig. 0. 

To obtain, for a given number of quarks exchanged, the total contribution of each 
graph, we must multiply the matrix element in eq.([TT|) by the number of ways of 
choosing the interacting pair without changing the structure of the graph and by the 
corresponding weight of the quark exchange operators in eq. @j. After multiplying 
these factors, the matrix elements in SFC space can be written as 

A + Bix + C&. (12) 

In Table |2], the values A, B and C for each type of graph and each number of quarks 
exchanged are shown together with the multiplying factors. 

It is known that the OGEP is not sufficient to describe the medium-range part 
of the baryon-baryon interaction, and that an attractive contribution arising from 
the meson cloud is also necessary. An effective meson exchange potential (EMEP) 
is introduced here as a force between H-dibaryons in order to simulate this meson 
contribution, although there are several works which introduce the meson exchange 



interactions as a force between quarks |T2|, 18|. The EMEP has previously been used 
to obtain a realistic description of the baryon-baryon interaction ||10|| . We extend the 
EMEP to describe the flavor-independent scalar meson exchange interaction between 
H-dibaryons. There are more elaborate meson exchange potentials |T2[, where other 
scalar and pseudoscalar meson exchanges are included, in addition to the simple flavour- 
independent scalar meson exchange used here. We shall argue later that the neglect of 
other meson exchanges will not substantially affect the qualitative conclusions obtained 
in this paper. 

Here we introduce an EMEP which is assumed to be flavour-independent. Since the 
EMEP is expressed as a function of the relative coordinate r between the H-dibaryons, 
we here explain how we add the EMEP to the RGM equation in fllO|) . First, we rewrite 
the RGM equation so that it has the same form as the usual Schrodinger equation, 
albeit with a non-local Hamiltonian. After introducing the renormalized H H relative 
wave function, 

X(r) = Jdr'N^ M (r,r') X (r'), (13) 



and the renormalized Hamiltonian, 



H(r,r>) = J dr"dr"'N^(r,r")[H RGM (r",r"') - E iut N RGM (r",r'")}N R ^(r'",r% 



the RGM equation is equivalent to the following non-local Schrddinger equation: 



Here E- mt is the internal energy, and E cm is the energy of the relative motion in the 
center-of-mass system. Then we add the EMEP, V(r), to the renormalized Hamiltonian 
kernel H(r,r'), i.e., H(r,r') — > H(r,r') + V(r)5(r — r'). For the EMEP, we take a 
gaussian-type function for simplicity, 



The parameters of V(r), ohh and Vohh; are obtained as follows. In the same manner 
as in ref. ||, we employ the direct convolution of the EMEP determined from the 
scattering length and effective range for the 1 Sq state of the nucleon-nucleon interaction, 
and obtain a H H = 1-00 fmQ The parameter Vohh is fixed so that the possible range 
of the H-dibaryon mass is consistent with the double hypernucleus event found at 
KEK 0, which is interpreted as A^Be or \\B with a AA binding energy of 8.5 ± 
0.7 or 27.6 ± 0.7 MeV, respectively. The corresponding range obtained for Vohh is 
- 1 096 (AA threshold) > Vohh > -1227 MeV ( A 3 A B) [-1136 MeV (i° A Be)]. On the 
other hand, Vohh = —1302 MeV is obtained by the direct convolution of the EMEP 
determined from the nucleon-nucleon interaction, which gives a rather large binding 
energy of the H-dibaryon from the AA threshold, 43.4 MeV. The values obtained for 
Vohh are summarized in Table [| 

3 Numerical results and properties of the H— H in- 
teraction 

In this section, we show the numerical results obtained by solving the RGM equation 
including the EMEP with the parameters explained in the previous section. 

1 The corresponding effective scalar meson mass [i c g can be obtained by equating the Fourier 
transform of the gaussian potential with that of the Yukawa-type for small momentum transfer, i.e., 
H e f{ = 2/anH — 400 MeV, which is not very far from the empirical cr-meson mass. 



(14) 




(15) 



V(r) = VoHHexp(-r 2 /a H H)- 



(16) 
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In Fig|| the equivalent local potentials of bound and scattering states and the adia- 
batic potential for Vohh = —1227 MeV are shown in order to appreciate the qualitative 
features of the H-H interaction, first. 

The equivalent local potential Veqlp(^) 13 i s defined as a potential which satisfies 
the following radial Schrodinger equation: 

h 2 d 2 _ 

- ^TT^oM + VmM r )Xo{r) = E cm Xo{r), (17) 

where \x is the reduced mass and Xo( r ) is the renormalized wave function obtained by 
solving the RGM equation with the non-local potential. The local potential obtained 
by this procedure is called the trivially equivalent local potential. We will discuss 
another way of defining the equivalent local potential later. 

The adiabatic potential 14a (r) is often used in order to get an idea in which regions 
of space repulsion and where attraction dominates. It is shown by the solid line in Fig|| 
The adiabatic potential is the ratio of the diagonal element of the GCM (Generator 



Coordinate Method: see refs. [10], |TT|) Hamiltonian kernel to the GCM normalization 
kernel. 

Kd(r) = ^ GCM K } - 5 GCM i r,r ! • (18) 



Nqcm (?",?"■) iV GC M (r,r) 
Here the value at infinite H-H separation is subtracted, so that the I4a(r) goes to zero 
as r —>■ oo. The GCM normalization and Hamiltonian kernels are the overlap and the 
expectation value of the Hamiltonian H with a (Os) 6 configuration for each H-dibaryon, 
where the origins of the two H-dibaryons are separated by a distance r. Because the 
origin is not the center-of-mass of the H-dibaryon, the adiabatic potential represents 
a smeared H-H potential, where the center-of-mass of each H-dibaryon fluctuates. 
Although the local and the non-local parts of the potential are smeared out, the range 
of the fluctuations is sufficiently small to reflect the short-range repulsion due to the 
effect of the exclusion principle and the CMI. 

We see from Fig.§ that the CMI leads to a strong short-range repulsion. The short 
range part of the EQLP, which is due to the CMI has an energy-dependent repulsion, 
i.e., the second bound and scattering states feel a stronger repulsion at short distances 
than the ground state. This energy-dependence is the typical one also seen in the 
nucleon-nucleon potential ||13|| . This behavior is shown to be partly due to the non- 
locality of the potential fl4j| . In the region below r < 2 fm, the contribution from the 
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OGEP has not yet died out and the EQLP still differs from the EMEP in the medium- 
range although the EMEP dominates there. Furthermore, the energy dependence seems 
to be reversed compared to the one at short distances. Namely, at short distances, the 
low energy state has a weaker repulsion than the high energy state, whereas at long 
range, the repulsion for the low energy state becomes stronger than the one for the high 
energy state. If we employ the usual procedure |H| to obtain the energy- dependent 



local potential from the non-local one, the ordering of the energy-dependence would 
be the same for short distances and long distances. 

Our results for the energy-dependence of the H-H interaction can be understood in 
the following way. In the interacting region, the wave function Xo( r ) for the attractive 
non-local potential is suppressed compared with the wave function Xedl(^) determined 
from the energy dependent local (EDL) potential. 



~ / \ m*{r) _ 

Xo{r) = \ Xedi», 

V m 

where m*(r) is the r-dependent effective mass J[5]. This is known as the Perey-Buck 
effect [|l6j. In order to obtain the energy-dependent equivalent local potential, we 
should use Xedl(^)- In the actual calculation, however, we have employed Xo( r ) instead 
of Xedl(^)- Then the kinetic energy given by the second derivative of the wave function 
Xo(?") differs from the one for XediX 7 ")- This is the main reason for the different energy 
dependences at short and long distances. 

The binding energy of the H-H system as a function of the parameter Vohh is shown 
in Fig.[|. The H-H system is a strongly self-binding system in the range, —1096 > 
Vohh > —1227 MeV, and its binding energy in the ground state is 100~170 MeV. Even 
a second bound state, which is also an S-wave state, appears in the range Vohh ^ —1200 
MeV. This is also seen from the divergent behavior of the H-H scattering length a as a 
function of Vohh shown in Fig.|5], and from the change of the S-wave phase shifts with 
Vohh shown in Fig. |6|. 

Fig.^ shows the equivalent hard-core radius ah, which is defined as 

a h = ~ (19) 

where 5 is the phase shift, and E cm = h 2 k 2 /2fi. When k goes to zero, the value ah 
coincides with the scattering length. For comparison, we show ah for the case without 
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EMEP. For k > 2 fm -1 , ah is determined mainly from the repulsive core, and is of the 
order of 1 fm. 

The wave function of the lowest bound state, Xo(r)/r, is shown for V hh = —1227 
MeV in Fig.|| Because of the strong repulsive core, the peak of the wave function 
is around 0.7 fm in spite of the large binding energy of the H-H system. The wave 
function is clearly localized around the minimum of the EQLP. In order to see the 
effect of the renormalization, the unrenormalized wave function, Xo( r )/ r ; which does 

1/2 

not contain the factor iV RGM , is also shown. The differences between the renormalized 
and unrenormalized wave functions can be understood as follows. The square root of 
the normalization kernel is expanded in terms of the harmonic oscillator wave functions, 



Unimi with the size parameter y'2/36 as 

^RGM^ O = E VJ^Klm( r )Unlm(r')- (20) 



nlm 

Here, /^tv is the eigenvalue of the normalization kernel, 

/2\ N 49 41 

"^-Hs) + i6b) -a**™- (21) 

where N = 2n + I. The angular momentum I is zero for the S-wave interaction. If 
ji n = 1 for all N in the expansion of eq.(|(J), the left-hand side becomes ^-function, 
5(r — r'). The deviation of fx^ from 1 is due to the quark exchanges between the 
H-dibaryons. The eigenvalues [in are || for n = and ^ for n = 1, respectively. The 
large contributions from n = and n = 1 cause the enhancement of the amplitude 
around the origin and the suppression at intermediate range by the renormalization. 

These results are obtained by employing the simple flavor-singlet EMEP. Therefore, 
further examinations are required especially concerning the validity of the flavor inde- 
pendent EMEP. Here, we wish to make a brief comment on the use of the EMEP which 
is taken to simulate a flavor singlet scalar meson exchange. Recently, more elaborate 



nucleon-hyperon interactions based on the SU(6) quark model [12] have been suggested. 
In ref. [|T2|, flavor-octet scalar meson exchanges are also included, in addition to the 
flavor-singlet scalar meson. Although we do not take into account the octet mesons in 
the present calculation, we qualitatively discuss their effect on the H-H system. The 
potential due to the flavor-octet scalar-meson exchange has the following structure: 

V(r) = f 1 ■ f 2 V(r), 
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where f 1 and f 2 are flavor SU(3) operators for baryons 1 and 2, respectively. This form 
is based on the SU(6) quark model. In general, there are two types of meson-baryon 
couplings, which are called F and D. It is known that the F type coupling, which gives 
the above flavor- dependence for the baryon-baryon interaction, is dominant for scalar 
meson exchanges |12|| . The expectation value of the flavor operator for the H-dibaryon 
is given as follows. 

<H\f l - f 2 | H >= + f 2 f -fl- f 2 } = -12. 

Here, we have used that < (/1 + /2) 2 >= for the H-dibaryon and that the H-dibaryon 
consists of flavor-octet baryons for which < f \ > is 12. Thus, the contributions from the 
flavor-octet scalar-meson exchanges in the H-dibaryon are repulsive because the radial 
part V(r) is always attractive for scalar meson exchanges. As a result, if we include the 
flavor-octet mesons and fix the strength of the EMEP from the binding energy of the 
H-dibaryon, the strength of the flavor- singlet part becomes more attractive than in the 
present calculation without the flavor-octet mesons. For the H H system, the flavor- 
octet meson exchanges do not contribute because the H-dibaryon is a flavor SU(3) 
singlet. Therefore, if the flavor-octet mesons were included, the H-H system would be 
more strongly bound. 

For the nucleon-nucleon system, octet scalar meson exchange is attractive in the 1 S'o 
channel. We expect, therefore, that the strength of the flavor-singlet part of the EMEP 
derived from the nucleon-nucleon potential is decreased compared to the case without 
octet mesons. On the other hand, as discussed before, the strength of the flavor- singlet 
part of the EMEP as derived from the binding energy of the H-dibaryon is increased 
when the octet mesons are taken into account. Thus, the strengths derived from 
the nucleon-nucleon scattering data on one hand, and from the H-dibaryon binding 
energy on the other hand, would become similar when we include the octet scalar 
meson exchanges. Recall that without the octet scalar mesons, the EMEP derived 
from nucleon-nucleon scattering was more attractive (—1302 MeV) than the EMEP 
derived from the limits on the binding energy of the H-particle (—1227 MeV ~ —1096 
MeV). 

Finally, we comment on the effect of pseudoscalar meson exchange. The pseu- 
doscalar meson exchanges do not contribute to the H-H system. The contribution 
from pseudoscalar meson exchanges to the binding energy of the H-dibaryon is known 
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to be about 5 MeV |T7|, |18|]. Therefore, even if we take into account the pseudoscalar 
meson exchanges, there is no drastic change for the binding energy of the H-H sys- 
tem. Thus, the conclusion that the H-H system has a strong medium-range attraction 
leading to a bound state remains unchanged. 

A bound state of two H-dibaryons is also obtained in the Skyrme model[|T^]. The 
binding energy of the "tetralambda" state is E-q = 15 ~ 20 MeV, which is smaller than 
the value obtained in our calculation. 

Finally, we comment on the implications of our results for the occurrence of H- 
matter. Tamagaki || discussed the possibility of a phase transition from neutron 
matter to H-matter at a density which is 6~9 times greater than the normal nuclear 
density po- in Tamagaki's calculation, the H-dibaryons interact via a hard core po- 
tential and an attractive square well potential outside the core. In their work, the 
depth of the attractive potential was assumed to be so weak that it can be treated as a 
perturbation. The depth of the square well potential has been determined as the value 
where the scattering length becomes zero for the first time, when the strength of the 
attractive potential is gradually increased, i.e. as the limiting strength that still gives 
a positive scattering length. The pertinent strength parameter corresponding to our 
EMEP can be obtained from Fig. [|. It is found to be Vohh = —638 MeV. Thus the 
attraction used in Ref. || is much weaker than the one used in the present calcula- 
tion. If the attractive H-H potential is indeed as strong as the EMEP employed in the 
present calculation, the critical transition density beyond which H-matter formation is 
energetically favorable may be appreciably lower. However, a quantitative estimate is 
not simple because of the inapplicability of perturbation theory under a strong attrac- 
tive potential as in the present calculation. A recalculation of the critical transition 
density for H-matter formation is beyond the scope of the present work (see however 
ref. [EDI). 



4 Summary 

We have investigated the interaction between two H-dibaryons in the quark cluster 
model with a one-gluon-exchange and quadratic confinement potential between con- 
stituent quarks. In addition, a phenomenological attractive meson exchange potential 
between H-dibaryons is used. The parameters employed for the calculation are de- 
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termined so as to reproduce the octet and decuplet baryon ground state masses and 
from the requirement that the H-dibaryon mass is consistent with the observed binding 
energy of two A's in a double hypernucleus. These parameters were also used in our 
previous work on the N-H interaction H ||. In our calculation, the effective scalar 
meson exchange potential is assumed to be a flavor SU(3) singlet. We have presented 
a qualitative discussion on the effects of octet scalar mesons and pseudoscalar mesons. 
We conclude that the following properties of the H-H system remain unchanged even 
if we include the octet scalar and pseudoscalar mesons. 

The main properties of the interaction between H-dibaryons can be characterized as 
a short-range repulsion due to Pauli blocking and the color-magnetic interaction, and a 
medium-range attraction due to flavor singlet scalar meson exchange. As a result of the 
present RGM calculation, we have obtained a strongly-bound state of two H-dibaryons 
with a separation about 0.8 ~ 0.9 fm due to the strong repulsive core. 

The present results suggest that the critical density at which the formation of H- 
matter becomes energetically favorable may be lower than the value p cr j t = 6po found 
in Ref. ||. This would have interesting astrophysical consequences. In the framework 
of the Walecka model, it has recently been shown [^0| that with the present H-H 
interaction, H-matter is unstable against compression. Thus, if the central density 
of a massive neutron star exceeds the critical density for H-matter formation, the 
energetically favorable compression of H-matter provides a possible scenario for the 
conversion of a neutron star into a strange quark star. 
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Figure 1: Quark exchange diagrams between two H-dibaryons 
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Figure 2: Ten types of two-body interaction with quark exchange 
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Figure 3: The equivalent local potential (EQLP) for bound states and a scattering 
state of the H H system, and the adiabatic potential for V oim = —1227 MeV. The 
solid curve denotes the adiabatic potential, the dashed curve denotes the EQLP of the 
lowest bound state, the chain curve denotes the EQLP of the second bound state, the 
double-dotted chain curve denotes the EQLP of the scattering state with center-of-mass 
energy E cm = 50 MeV. 
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-1400 -1200 -1000 -800 -600 

Vohh (MeV) 

Figure 4: The binding energy of the H-H system as a function of the strength parameter 
of the effective meson exchange potential Vohh- Both bound states are S-wave states. 
The vertical dotted lines show the range of Vohh from —1227 to —1096 MeV. 
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Figure 5: The H-H scattering length a as a function of the strength parameter of 
the effective meson exchange potential Vohh- o — —1-26, —4.91, 6.12 and 3.51 fm for 
\/ 0HH = -1096, -1136, -1227 and -1302 MeV, respectively. The vertical dotted lines 
show the range of Vohh from -1227 to -1096 MeV. 
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Figure 6: The HH S-wave phase shift 5 as a function of the relative wave number k 
for several values of the strength parameter of the effective meson exchange potential 
Vqhh- 
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Figure 7: The equivalent hard core radius as a function of the relative wave number k 
for several values of the strength parameter of the effective meson exchange potential 
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Figure 8: The wave function of the lowest S-wave bound state for V hh = —1227 MeV. 
The solid curve denotes the wave function without the renormalization Xo( r )/ r and the 
dashed curve denotes the renormalized wave function Xo( r )/ r - The root-mean-square 
radii of the H-H system are 0.92, 0.89, 0.84 and 0.81 fm for \/ 0H h = -1096, -1136, 
-1227 and -1302 MeV, respectively. 
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Table 1: Quark model parameters 



b (fm) 


m (MeV/c 2 ) 


a s a c (MeV/fm 2 ) 


6 


6 


0.6 


300 


1.394 33.0 


0.6 


0.1 



Table 2: The matrix elements of spin-flavor-color parts 



number of exchanged quarks 


overlap 


interaction type 




• (TjXi ■ 


A, 




Xi ■ Xj 




(factor) 






(factor) 


A 


B 


C 


A 


B 


c 






1 


(15) 


K 
o 


99 


_Q 
o 


_7 
i 


22 
3 


5 
3 






2 


(15) 


5 


22 


-3 


-7 


22 

3 


5 
3 






3 


(0) 
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(0) 


o 


o 
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o 


o 
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(1) 
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o 
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o 


o 






10 


( 0) 
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o 


o 


o 


o 


o 






1 


(101 


20 
3 


88 
3 


4 


28 
3 


88 
9 


20 
9 






2 




20 
3 


88 
3 


4 


28 
3 


88 
9 


20 
9 








(5) 


10 


44 




14 


44 


10 






3 


3 


3 


2 


3 


9 


9 


1 
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(5) 


10 

3 


44 
3 


2 


14 
3 


44 
9 


10 
9 




-2 


5 


(5) 


10 

3 


44 
3 


2 


14 
3 


44 
9 


10 

9 


you) 




6 


(5) 


10 

3 


44 
3 


2 


14 
3 


44 
9 


10 
9 






7 


(25) 


89 


11 


139 


683 


341 


19 






16 


24 


48 


144 


72 


16 
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(1) 


64 
3 
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32 
3 


64 
9 


o 


32 
9 
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(0) 


o 
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o 


o 


o 






10 


(0) 


o 


o 


o 


o 


o 


o 






1 


(6) 


13073 


2335 


3971 


11851 


21623 


6965 






1296 


72 


1296 


1296 


1944 


3888 






2 


(6) 


13073 


2335 


3971 


11851 


21623 


6965 






1296 


72 


1296 


1296 


1944 


3888 






3 


(8) 


7 


253 


137 


233 


253 


223 








8 


12 


24 


24 


36 


72 


2 




4 
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7 


253 


137 


233 


253 


223 
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12 
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24 


36 


72 




49 
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7 


253 
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233 


253 


223 




16 
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12 


24 


24 


36 


72 


(225) 
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8 


253 
12 


137 
24 
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24 


253 
36 


223 
72 






7 
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36281 


6593 


21437 


18347 


7357 


5495 






3888 


1944 


3888 


1944 


972 


1944 






8 


(4) 


8227 


1786 


3667 


23611 


2093 


10759 








243 


243 


243 


1944 


972 


1944 






9 


(1) 


3953 


499 


269 


1675 


4187 


419 






648 


36 


648 


648 


972 


1944 
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6545 


853 
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2527 


677 
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864 
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1728 
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87 
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32 
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32 
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32 
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3 

32 
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16 
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32 
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32 


57 
16 


43 
32 






7 


(9) 


33185 


5239 


13973 


13611 


863 


1837 






2592 


1296 


2592 


2592 


432 


864 






8 


(9) 


31319 


5845 


12737 


8031 


23 


1327 








1296 


648 


1296 


1296 


216 


432 






9 


(3) 


2131 


8575 


1091 


3895 


2513 


691 






576 


864 


1728 


1728 


864 


1728 






10 
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2131 


8575 


1091 


3895 


2513 


691 






576 


864 
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1728 


864 
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Table 3: Strength parameter of the effective meson exchange potential 



Vohh (MeV) 2M A - M H (MeV) 

-1096 ~ AA threshold 

-1136 8.5 A ° A Be 

-1227 27.6 A 3 A B 

-1302 43.4 
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